Updating Boer-Mulders functions from unpolarized pd and pp Drell-Yan data 
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We extract the Boer-Mulders functions for the proton by combining the unpolarized pd and pp Drell-Yan 
data measured by the E866/NuSea Collaboration by the assumption that the cos 20 asymmetry is from the 
Boer-Mulders function. Using the extracted Boer-Mulders functions, we present the predictions for the cos 20 
asymmetries in future pp experiments at J-PARC and pp experiments at PANDA and PAX. 
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I. INTRODUCTION 

In a recent work iH] the Boer-Mulders functions for the 
proton have been extracted from the cos 20 angular asym- 
metry measured in the unpolarized p d Drell-Yan process by 
the E866/NuSea Collaboration lH at FNAL. Much earlier, 
the first cos 20 asymmetries of dilepton production had been 
measured two decades ago by the NAIO [3J Collaboration 
and E165 |4] Collaboration, but for the 7r-nucleus Drell-Yan 
processes, showing that the magnitude of the asymmetries is 
around 30% at most. This substantial cos 20 angular depen- 
dence that violates the so-called Lam-Tung relation js]] pre- 
dicted by perturbative QCD, belongs to the remaining chal- 
lenges which need to be understood from QCD dynamics. 
Several attempts have been made to interpret these data, in- 
cluding QCD vacuum effects |0, ^ and higher-twist mech- 
anisms (HI 1^. Furthermore, in the last decade significant 
efforts have been put forward on the understanding of the 
cos 20 angular dependence from the view point of the trans- 
verse momentum dependent (TMD) Boer-Mulders function 
h^{x,p^) It was shown [TT] by Boer that the angular 

dependence can be related to the product of two functions 
hf, each of which comes from one of the incident hadrons. 
The Boer-Mulders function describes a correlation between 
the transverse spin and the transverse momentum of a quark 
inside an unpolarized hadron. Despite the naive T-odd nature 
of the Boer-Mulders functions and of their chiral-even part- 
ner the Sivers functions [ 12], it has been shown that they can 
originate from inital/final state interactions ifTsi [fill between 
the struck quark and the spectator of the nucleon, which are 
important to ensure the gauge invariance of the TMD distri- 
bution functions iflsl - flTll . Recently a lot of theoretical stud- 
ies and phenomenological analysis iUsI - IstIi on Boer-Mulders 
functions have been carried out. 

The first measurement of the cos 20 asymmetry in 
the nucleon-nucleon interacting Drell-Yan process by the 
E866/NuSea Collaboration makes the attempt on extracting 
the proton Boer-Mulders function possible. In Ref. Hi] we 
parameterized the Boer-Mulders functions for m, d, u and d 
quarks, which were fitted to the E866/NuSea p d data, based 
on the assumption that the cos 20 asymmetry comes only from 
the Boer-Mulders effect in the region <K Q^, where qj and 
Q are the transverse momentum and the invariant mass of the 
lepton pair In our parametrization the transverse momentum 



dependence of hf{x,pl.) was modeled by the Gaussian ansatz: 
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The parameter pl^ describes the Gaussian width of the trans- 
verse momentum distribution. The x-dependence of /i|^(jc) 
was further parameterized to be proportional to the unpolar- 
ized parton distribution h^'^(x) - x"^ (1 - x) f^'{x), where the 
parameter c is the same for all flavors. The large-x behavior 
of h^'^ix) compared with f^{x) was taken into account by the 



factor( \ - x) following the argument in Ref. 13811 . 

More recently the E866/NuSea Collaboration reported new 
measurements [39] on the cos 20 asymmetry in the unpolar- 
ized pp Drell-Yan process. The overall magnitude of the 
cos 20 dependence for pp processes is consistent with, but 
slightly larger than that of pd processes. The new pp data, 
besides the previous p d data, will provide further information 
and constraint on the shape of the Boer-Mulders functions for 
diff'erent flavors. Therefore there is the need to perform a new 
extraction of Boer-Mulders function in the presence of new 
pp data. In this work we combine the previous pd data and 
the new pp data to extract Boer-Mulders functions for the pro- 
ton to update our previous results. Furthermore, in the new fit 
we include x/-dependent and ^-dependent data that have not 
been applied in the previous extraction. We then apply our ex- 
tracted Boer-Mulders functions to predict the cos 20 asymme- 
tries in future pp experiments at J-PARC and pp experiments 
at PANDA and PAX. 



II. DESCRIPTION OF cos 20 ASYMMETRIES IN TERMS 
OF BOER-MULDERS FUNCTIONS 



The angular differential cross section for unpolarized Drell- 
Yan processes has the general form 
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where 9 and are, respectively, the polar angle and the az- 
imuthal angle of dileptons in the CoUins-Soper frame ll40ll . 
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The coefficients A, fi and v do not depend on ttiese angles, and 
for scattering that has azimuthal symmetry their values are 
ju = V = 0. 

This angular distribution has been measured in muon pair 
production by pion-nucleon collisions: n^N — > fi^fi^X, with 
N denoting a nucleon in deuterium or tungsten, and for a 
TT" beam with energies of 140, 194, 286 GeV B] and 252 
GeV [4]. The experimental data show large values of v, 
near 30%. The most recent measurements of the angular dis- 
tribution were performed by the E866 Collaboration |2], in 
pd Drell-Yan processes at 800 GeV/c. The measured v is 
about several percent, a result which can not be explained by 
leading-twist collinear factorization ' in QCD. As proposed 
by Boer ifTTIl . the non-zero cos 20 term can be produced by 
the product of two hfs, each coming from one of the two in- 
cident hadrons. 

The leading order unpolarized Drell-Yan cross section ex- 
pressed in the Collins-Soper frame ll40ll is ifTHl 
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where = and qr are the invariance mass square and 
the transverse momentum of the lepton pair. The vector h - 
qr/Qr- We have used the notation 



r[ff] = j d'pTd'krd'ipT + kT- qT)f{xl,p'T).f(X2, k'j).(A) 

The first term in Eq. ([3]) is azimuthal independent, while the 
second term has a cos 20 azimuthal dependent term which 
contributes to the asymmetry v. 

In the case of the pd and pp Drell-Yan processes, the 
cos 2(j) asymmetry can be expressed as 



VpdiXy,X2,qT) 
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where 



(7) 



Fpd(xi,X2,qT) = 2r[i2h-pTh-kT-pT-kT)(elhf" 

+ ejhf'Xhf + hf')] + (q^q), 
Gpdixi , X2, qr) = T [{elf- + elff)(fl' + ff)] + (q ^ q), (8) 
FppixuX2,qT) = 2r[i2h-pTh-kT-pT-kT)ielhf"hf 



' In has been shown in Ref. (aT] that the cos 20 asymmetries can be ex- 
plained by taking into account the twist-three quark-gluon correlations in 
collinear factorization which is consistent with the Boer-Mulders effect in 
the TMD factorization approach. 
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FIG. 1: The ^7-dependent cos 20 asymmetries for unpolarized pd 
(left) and pp (right) Drell-Yan process calculated from our fitted re- 
sults. The region between the two dotted line corresponds to the 
uncertainty of the parameters. Data are from the FNAL E866/NuSea 
experiments. Here we exclude data with pr > 2 GeV. 



+ elhf^hf')] + iq^q), (9) 
Gpp(xuX2,qT) = r[(elf''fl'+elffff)]+(q^q). (10) 

For distribution functions for deuteron, we have used the 
isospin relation: 



u I deuteron 
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In Ref. fit] we have parameterized the transverse momentum 
dependence of Boer-Mulders functions with a Gaussian form 
as follows 
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The X dependence for m, d, u and d quarks is parameterized, 
as follows 



h^^x) = //,x^(l-x)/jV). 



(13) 



The above parametrizations, with 6 parameters, have been ap- 
plied Xofilpd Drell-Yan data measured by E866/NuSea Col- 
laboration. In the fit the f j--dependent and xi/2-dependent 
cos 20 asymmetry data were used. The fitted result was em- 
ployed to predict the x/^-dependent and Q-dependent cos 20 
asymmetries which were compared with the coiTesponding 
data. 

Recently the E866/NuSea Collaboration reports measure- 
ments 1391 of the cos 20 asymmetries on unpolarized pp 
Drell-Yan processes at Ep - 800 GeV. The new pp data, to- 
gether with the previous p d data, will provide further infor- 
mation on the shape of the Boer-Mulders functions for differ- 
ent flavors. In this paper we will combine the previous p d data 
and the new p p data in the fit. Further more we will include 
xy^-dependent and Q-dependent data in our fit. To do this we 
will parameterize the Boer-Mulders functions as in Eqs. (fT2b 
and ( fT3] l. but changing the form slightly. 

In our previous fit we modeled the jc-dependent behavior 
of h\''{x,p\) at small x as x"^ compared with f^ix), and we 
assumed the value of c to be flavor independent, as shown in 
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FIG. 2: The X| -dependent cos 20 asymmetries for unpolarized pd 
(left) and pp (right) Drell-Yan process calculated from our fitted re- 
sults. The region between the two dotted line corresponds to the 
uncertainty of the parameters. Data are from the FNAL E866/NuSea 
collaboration. 




FIG. 4: The -dependent cos 20 asymmetries for unpolarized p d 
(left) and pp (right) Drell-Yan process calculated from our fitted re- 
sults. The region between the two dotted line corresponds to the 
uncertainty of the parameters. Data are from the FNAL E866/NuSea 
collaboration. 
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FIG. 3: The X2-dependent cos 20 asymmetries for unpolarized pd 
(left) and pp (right) Drell-Yan process calculated from our fitted re- 
sults. The region between the two dotted line corresponds to the 
uncertainty of the parameters. Data are from the FNAL E866/NuSea 
collaboration. 
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FIG. 5: The Q-dependent cos 20 asymmetries for unpolarized pd 
(left) and pp (right) Drell-Yan process calculated from our fitted re- 
sults. The region between the two dotted line corresponds to the 
uncertainty of the parameters. Data are from the FNAL E866/NuSea 
collaboration. 



( fTSl l. Now with more data available, we are able to release 
this constraint to replace c as Cq, depending on flavor. Sec- 
ondly we model the large jc-dependence of the Boer-Mulders 
functions by (1 - x)*, different from our previous fit in which 
the large x dependence is 1 - x. Therefore we have the new 
parametrizations for h^'' (q - u,d,u and d) as follows: 



(14) 



The first p^-moment of Boer-Mulders function is defined 



as 



h^'^'^^ix)^ j(fpT^hf(x,pl) (15) 

From Eqs. (fTSl i and ( fT4l l one can calculate from our 

parametrization as 
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With the Gaussian form for the pj dependence of Boer- 
Mulders functions and the unpolarized TMD distribution 
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(17) 



the transverse momentum integrations in Eqs. (|7]l - ( fTOl i can 
be deconvoluted and the results are: 



Fpd{xuX2,qT) = Fpd(xuX2) 



Fpp(xi,X2,qT) = Fpp(xi,X2) 



Gpp(xi,X2,qT) - Gpp(xi,X2) 

where the xi/2 dependent parts are 

Fpd(xuX2)^ {4hf"(xi)+ hf'(xi)) 
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X {hf"(x2) + hf\x2)) + (q ^ q) 

= (l-Xif(l- X2)\^Hi x\" f^(xi) /f(x2) 

+ H2X\' f({xi)x'^ff(x2) +4H,x\"f{'(xi)xl'ff(x2) 

+ (H1H2/H3) x\' ff(xi)xl' /ffe)) + (q 



1 ji\-^ij-^2 ./iv-^z// ' vi ' q)^ (22) 

FppixuX2) ^4hf" (Xi)hf"(x2) + hf''(xi)hf'(x2) + (q^q) 
b 



= (1 - xi)" (1 - X2fiAH, x\' /i"(x,) /i"(x2) 



+ H2X\' f,(xi)x^ f^(x2)) + {q ^ 
X (/f (X2) +/f(X2)) + (? ^ 

Gp„(x,,X2) = 4/i"(x,)/f(x2) + /f(xi)/f(x2) + ^ ^) 



where Hi = //„//,-„ //2 = //rf//,?, H3 = //„//j and 
H\H2IHt, - HjjHd. Since //,„ Z/^/, //,-, and //s always appear 
as products of two of them, we will apply Hi , H2 and as 
the parameters in the fit. Therefore the actual number of free 
parameters is reduced to 9. 

The qr-, xi - and X2-dependent cos 2<p asymmetries in unpo- 
larized p d and pp Drell-Yan processes can then be expressed 
as 
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where the subscript A^A^ denotes p d and pp. 

One can also express the cross-section of the Drell-Yan pro- 
cess, depending on Feynman xf and the mass of the lepton 
pair Q as 
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with 
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Therefore the Xf- and Q-dependent cos!] 
be expressed as 



xi X2 s. (30) 
b asymmetries can 
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FIG. 6: The first p^-moments of Boer-Mulders functions for u, d, u 
and d quarks for = I GeV^ by solid lines, the shadows depict the 
variation ranges of xh'^^^'"'{x) allowed by the positivity bound. The 
dashed lines show x/j' (x). . 



III. FITTING BOER-MULDERS FUNCTIONS TO THE 
UNPOLARIZED E866/NUSEA AND pp DATA 

The E866/NuSea Collaboration measured the cos 2ip asym- 
metry ies Vpd and Vpd \s Qt, xi, X2, xf and irip^ in the following 
kinematical region: 

4.5 GeV <Q<9 GeV and 10.7 GeV < g < 15 GeV, 
^j-<4GeV, 0.15 < xi < 0.85, 0.02 < X2 < 0.24. 



In the following we apply the theoretical expressions ( 

), (|3TT i and ( [32] i to fit the unpolarized p d and pp Drell-Yan 

cos 2<p asymmetry data 101391]. The Boer-Mulders effect to the 
cos 2<p asymmetry is supposed to apply in the region where qj 
is not large. At large qr, the higher order perturbative QCD 
contributions ll42i l43ll might be important. Therefore we ex- 
clude the data with qj > 2 GeV in our fit. For the parton 
disti-ibution /f (x) we adopt the MSTW2008 LO set H. We 
choose the Gaussian width for f^jx, pi) as - 0.25 GeV^, 
following the value given in Refs. ll45ll46ll . The best fit results 
and the errors for the parameters are as follows: 



Hi 

Cu 



0.62_o29' 


H2 


- 1 45+1-30 


H3 


n /;i +0.50 


0-63_o.2i' 


Cd 


- 47+"-36 

- U.t/_o.39> 


Cu 


- 07+" "'' 

- U.U/_() 05, 


^■'-'-0.52' 


bo 


- u.i /_() 14, 


2 

Phm 


- 173+0 027 



(33) 



The of this fit is 35.95 for 52 data points, resulting 
X^/d.o.f = 0.84. In Figs.[T]|2l|3]|4]and|5]we show the qr-, xi-, 
X2-, Xf - and 2-dependent cos 20 asymmetries for unpolarized 
pd and pp Drell-Yan process calculated from our fitted re- 
sults and compare them with FNAL E866/NuSea data. The 
solid lines show the best fit results, and the regions between 
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FIG. 7: The cos 2(j) asymmetries vs Xf at J-PARC with proton beam 
energy Ep = 50 GeV. The region between the two dotted hne corre- 
sponds to the uncertainty of the parameters. 



the two dotted lines correspond to the uncertainty of the pa- 
rameter In Fig.[T]we also show the predictions for pj > 2GeV 
region from the Boer-Mulders effect. 

The possible range of coefficients Hg are obtained from 
the values of Hi, H2 and H^, by employing the positivity 
bound ^47*1 for h^'^(x,p^) for the entire x and pr regions: 



\pTh^''ix,pl)\ 
M 



<fi'(x,pl). 



(34) 



we assumed that it has the same behavior as that of the un- 
polarized distribution f'^. Recently there has been gro wing 
interest on performing next-to-leading order analysis 11481 - 
Is^] of single-spin asymmetries in Semi-inclusive deeply in- 
elastic scattering and Drell-Yan processes, especially on the 
Sivers asymmetry. Those studies can provide information on 
the evolution of the Sivers function. It is also interesting 
to study the evolution of the Boer-Mulders function. How- 
ever, to investigate the scale dependence behavior of the Boer- 
Mulders function and its impacts on the cos 2<p asymmetries 
is out of the scope of this paper. Secondly, we considered 
the Boer-Mulders effect as the dominant source for the cos 20 
asymmetries, which is reasonable in the region <K Q^. 
Since the data in E866/NuSea covers the kinematics regime 
4.5 GeV < Q < 9 GeV and 10.7 GeV < Q < 15 GeV, one 
expects that the Boer-Mulders effect dominates in the region 
qr < 2 GeV. At qr > 3 GeV, the higher order perturbative 
QCD contributions (4^, '43'] might be important. This can be 
seen from Fig.lT] which indicates that the predicted asymme- 
tries at large qr are small, while the size of the data in that 
region is substantial f^, '39']. Thirdly, the x^/d.o.f in the fit 
presented in this work is a little bigger than that in the previ- 
ous fit shown in Ref. 1 1]. This is because in Ref. ^ we only 
included 16 data points from the pd process, while now we 
have 52 data points in our new fit. And among those data we 
include Q- and .it/r-dependent data also. We have checked that 
the two data points in the large Xf region give rise to substan- 
tial contributions (20%) to the total in our best fit. 



We have 
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(35) 



The upper and lower limits for Hq are determined by the pos- 
itivity bound for h^''{x,p^). The central value for Hq shown 
above is obtained from the geometric mean values of the up- 
per and lower Umits for Hq-. H^^^ = ^H^^^H^^^. In our 
previous work 1 1 ] the variation range of Hq allowed by the 
positivity bound was described by the coefficient co, namely, 
that the substitution Hq — > u)Hq for q - u, d and Hq — > ^Hq 
for q - u,d will not change the result. In our new fit presented 
here, the range of o) is 0.48 < w < 2.1, and central values for 
Hq correspond to w = 1 . 

The positivity bound given in ( [34l ) implies 



(36) 



In Fig.|6]we show the first /7^-moments of the Boer-Mulders 
functions xh''^^^^''(x) for u, d, u and d quarks for - I 
GeV^ by solid lines, the shadows depict the variation ranges 
of xh^^^^''(x) allowed by the positivity bound. The dashed 

lines show ^^g^-^ (x ). 

Several comments on our fit are in order First, since we 
don't know the scale dependence of Boer-Mulders functions. 



IV. PREDICTIONS FOR FUTURE pp AND pp 
EXPERIMENTS 



We will then apply our extracted Boer-Mulders in previous 
section to predict the cos 2<p asymmetries in future pp and pp 
experiments. 

Drell-Yan process has been proposed at J-PARC llssll by pp 
scattering with proton beam energy Ej, = 50 GeV. with this 
lower beam energy than that at E866/NuSea, the measurement 
of cos 2(p asymmetry at J-PARC will provide complementary 
information on Boer-Mulders functions in a different kinemat- 
ical region. We estimate the Xf - dependent cos 2<p asymme- 
try at J-PARC by imposing the cuts < qr < I GeV and 
4 < g < 5 GeV from the fitted results in Eq. (I33]l directly, 
as shown in Fig.|7] The solid curve shows the result from the 
best fitted values, while the region between the dotted lines 
correspond to the parameters uncertainty. 

The GSI-PANDA experiment Q will run the Drell-Yan 
process with an unpolarized antiproton beam colliding with 
an unpolarized proton target at s - 30 GeV^. The PAX [55| 
experiment might also perform unpolarized pp Drell-Yan pro- 
cess at the energy 5 = 45 GeV- in the fixed target mode. The 
kinematical cuts we apply in the estimation of the asymmetry 
for those experiments are 2 < 2 < 3 GeV and < ^7- ^ 0.4 
GeV. To estimate the asymmetries in pp processes we need 
to use the Hq results given in (l35l l. rather than to the results 
for Hi (i - 1, 2, 3) in (|33] |. In Fig. |8] we show the predicted 
cos2(/> asymmetries v at PANDA for s - 30 GeV-, and PAX 
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FIG. 8: The cos 20 asymmetries vs Xf at PANDA for i = 30 and 
at PAX for i = 45 GeV^, respectively. The bands correspond to all 
possible values of a> in the allowed range 0.48 < co < 2.1. 



for 5 = 45 GeV^, as a function of Xf. The solid lines are cal- 
culated from the central value of Hg given in Eq. ( |35] ). and the 
bands correspond to all possible values of oj in the allowed 
range 0.48 < a> < 2.1. The asymmetry calculated from the 
central value of H^j in pp process is smaller than that in pp 
case. Therefore the measurement on the cos 2(f> asymmetry at 
PANDA and PAX might provide valuable examination on our 
extraction of Boer-Mulders functions. 



V. SUMMARY 

We have parameterized the Boer-Mulders functions for the 
proton by employing a Gaussian form for their transverse mo- 
mentum dependence. We then fitted our parameterizations to 
both the previous measured pd data and recent pp data on 
the unpolarized Drell-Yan cos 20 dilepton asymmetries from 
the E866/NuSea Collaboration. The basic assumption in the 
fit is that the cos 20 asymmetry is contributed by the prod- 
uct of two Boer-Mulders functions and that other contribu- 
tions such as perturbative QCD effects can be ignored in the 
region <K Q^. In our fit we included not only the qr-, 
xi- and jC2-dependent asymmetry data, but also Q- and Xf- 
dependent data. We applied our extracted Boer-Mulders to 
predict the cos 20 asymmetries in future pp experiment at J- 
PARC and pp experiment at PANDA and PAX. We found a 
smaller asymmetry for pp processes compared with that for 
pp processes, which might serve as a test on our extraction of 
Boer-Mulders functions. 
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